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Abstract 

In this paper, specific features of Sm magnetism in an intermetallic compound have been studied. 
For this purpose, a high quahty single crystal of SmPd2Al3 was grown and subjected to detailed 
measurements of specific heat, magnetization, AC susceptibility and electrical resistivity with re- 
spect to temperature and a magnetic field applied along the principal crystallographic directions. 
SmPd2Al3 magnetism was found to be strongly anisotropic with the easy-magnetization direction 
along the c-axis where the main magnetic features are concentrated. The a-axis response remains 
weak, paramagnetic-like, even in the magnetically ordered state. Ferromagnetism with Tc = 12.4 
K has been indicated by all the measured physical properties. At lower temperatures, three succes- 
sive order-order phase transitions have been observed on the temperature dependence of the specific 
heat as three anomalies: at 3.4, 3.9, 4.4 K, respectively. The low temperature magnetization data 
can be understood within a scenario that considers the antiferromagnetic ground state as being 
gradually destroyed through a series of four metamagnetic transitions at 0.03, 0.35, 0.5 and 0.75 
T, as detected in the 1.8-K magnetization data. The temperature dependence of the paramagnetic 
susceptibility below 200 K can be in the first approximation interpreted in terms of a Curie- Weiss 
law modified by temperature independent van Vleck contribution due to the low-lying first excited 
multiplet J = 7/2 being populated. At higher temperatures, involvement of the second excited 
multiplet J = 9/2 should also be considered. The experimental data are discussed together with 
the results of electronic-structure and crystal- field calculations from first principles, which were 
performed as an important part of the study for comprehension and explanation of the observed 
behavior of the SmPd2Al3 compound. 

PACS numbers: 75.30.-m; 75.40. Cx; 
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I. INTRODUCTION 



The localized -!^/-electron magnetism of lanthanide ions is by rule characterized by stable 
magnetic moments, which reflect the population of energy levels of the ground-state multi- 
plet determined by total angular momentum J and split by crystal fleld (CF) interaction. 
Paramagnetic susceptibility (x) well above the magnetic-ordering temperature obeys the 
Curie Weiss law. This behavior may be modifled at lower temperatures in the case of a 
considerable CF interaction. Sm represents an exception to this rule: Sm magnetism is not 
governed only by the ground state multiplet J = 5/2, since the flrst (and second) excited 
multiplet J = 7/2 (J = 9/2) has only an energy of 0.1293 eV (0.2779 eV)^ Therefore the 
temperature dependence of paramagnetic susceptibility at temperatures up to about 400 
K is influenced by populating the CF split levels of the flrst and second excited multiplets 
via the temperature-independent Van Vleck terms^, and does not obey the Curie- Weiss law. 
A typical 1/x vs. T plot for a Sm compound is strongly nonlinear and usually exhibits a 
maximum around 400 K, which is very sensitive to CF and exchange interaction^. Since the 
samarium magnetic moment at low temperatures is small {gJ = 0.71 /i^ for the Sm^"*" free 
ion) exchange interactions can be very peculiar and several magnetic phases may appear in 
the same material at different temperature and magnetic-fleld intervals, respectively. In the 
present paper we demonstrate some speciflc features of Sm magnetism from the results of 
our recent study of a high-quality single crystal of the compound SmPd2Al3. We focused 
our attention on measurements of magnetization (M), AC susceptibility (xac), speciflc heat 
(Cp) and electrical resistivity (p) as functions of temperature (T) and the external magnetic 
fleld (B). The anisotropy of magnetic properties has been determined from measurements in 
magnetic flelds applied along the principal crystallographic directions [100] and [001]. Along 
with our experimental work, we performed various theoretical calculations of the electronic 
structure and CF parameters and compared calculated results with experimental flndings. 
SmPd2Al3, as well as its i?£Pd2Al3 [RE = La, Ce, Nd, Gd) counterparts, adopt the same 
hexagonal PrNi2Al3-type crystal structure^. It is worth noting that the U isostructural ana- 
log, UPd2Al3, is a well-known heavy-fermion antiferromagnetic superconductor (7 = 150 
mJ K~^mole~^, Tc = 2 K, Tn = 14 K)^. The nonmagnetic analogue LaPd2Al3 becomes 
superconducting at 0.8 K^. The susceptibility of the Ce-, Pr-, Nd- compounds obeys the 
Curie- Weiss law at temperatures above 100 K. 
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Deviations from the Curie- Weiss law at temperatures below Tc are apparently due to 
CF splitting of the ground-state multiplet^, which is also reflected in the Schottky anomaly 
observed in specific-heat data collected for CePd2Al3 and PrPd2Al3^. CePd2Al3 exhibits 
heavy-fermion behavior at low temperatures (7 = 380 mJ K^^mole^^) with a Kondo effect;- 
nevertheless, it becomes antiferromagnetic below 2.8 KiS^^. PrPd2Al3 remains paramagnetic 
down to 1.5 K. Saturation of susceptibility below 15 K has been attributed to the CF effect 
yielding the ground-state singlet^. NdPd2Al3 becomes antiferromagnetic at Tn = 6.5 K as 
reported in^. The value of of NdPd2Al3 has been found to be varying in the interval from 
5.2 to 7.7 K for various samples, following a linear dependence on the lattice parameter a, not 
on c, which may be understood with the exchange interactions inside the ah plane dominating 
the magnetism in this compound^^. The most intriguing case in the i?iiPd2Al3 family seems 
to be SmPd2Al3. Results published so far on SmPd2Al3 have been obtained by measuring 
only polycrystals and point to a complex magnetic phase diagram with several magnetic 
phase transitions below 12 K-. The CF interaction in SmPd2Al3 can be reasonably estimated 
by considering the well-determined CF parameters for the NdPd2Al3 compound^*^ or by 
using CF calculations as obtained from first principles, as will be presented in this paper. 
Specific-heat measurements of SmPd2Al3 have been reported revealing anomalies at 3, 6 
and 12 K, indicating magnetic phase transitions*. From susceptibility data three magnetic 
phase transitions at 4, 4.3 and 12 K-, respectively, have been proposed whereas the electrical 
resistivity reflects strongly a magnetic phase transition at 12 K^, typical for a transition to 
ferromagnetic ordering, the peak- like anomaly around 3.5 K resembles a transition to a 
low-temperature antiferromagnetic phased. The magnetization curve measured at 2.5 K 
indicated two magnetic phase transitions (at approximately 0.3 and 2 T)^°. It is clear that 
SmPd2Al3 presents unique case of 4f magnetism and we attempt here to put this context 
together with previous results and Sm magnetism generally. 

II. EXPERIMENTAL AND COMPUTATIONAL DETAILS 

The single crystal of SmPd2Al3 was grown in a triarc furnace by the Czochralski pulling 
method from stoichiometric amounts of elements. The single crystal was a 20-mm long 
cylinder with a diameter of 3 - 4 mm. Crystal quality was checked by the Laue technique. 
A small part of the crystal was pulverized in an agate grinding mortar and X-ray powder 
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diffraction (XRPD) data were recorded on a Seifert powder diffractometer equipped with a 
monochromator to provide the Cu Ka radiation. The XRPD data were analyzed by means 
of the Rietveld profile procedure^^ using the program FullPro^. The final composition of 
the crystal was checked by EDX analysis on a FE-SEM Tescan. The Laue technique was 
also used for crystal orientation. Samples appropriate by shape for the magnetization and 
specific heat measurements were cut by a fine wire saw (South Bay Technology Inc, type 
810). The sample for the magnetization measurements had the form of a small beam (1x1 
X 1.5 mm^) with the rectangular planes oriented perpendicular to the crystallographic axes a 
and c, respectively. The specific-heat samples were small plates (1.5x1.5x0.2 mm'^) with the 
main-plate plane perpendicular to the a- and c-axis, respectively. The orientation of each 
prepared sample was subsequently revised again by Laue technique. The specific-heat, AC 
susceptibility and magnetization measurements were performed using the Quantum Design 
PPMS (Physical Property Measurement System) and the MPMS (Magnetic Property Mea- 
surement System), respectively. The specific heat was measured at temperatures from 1.8 to 
300 K in magnetic fields up to 9 T. The magnetization was measured in a temperature range 
from 1.8 to 400 K and in magnetic fields - 14 T. The AC susceptibility was measured from 
1.8 to 30 K in an AC magnetic field 0.03 mT with a frequency of 497 Hz. The magnetic field 
was applied along the a- and c-axis, respectively. To obtain direct information about the 
ground-state electronic structure and magnetic properties, we applied theoretical methods 
from first principles. The ground-state electronic structure was calculated on the basis of 
density functional theory (DFT) within the local spin density approximation (LSDA)^ and 
generalized gradient approximation (GGA)^. For this purpose, we used the full potential 
augmented plane wave plus local-orbitals method (APW -|- lo) as implemented in the latest 
version (WIEN2k) of the original WIEN code^^. The calculations were scalar relativistic 
and performed with the following parameters: non-overlapping atomic spheres (AS) radii 
of 2.8, 2.5 and 2.0 a.u. (1 a.u. = 52.91 pm) were taken for Sm, Pd and Al, respectively; 
the basis for the expansion of the valence electron states (less than 6 Ry below the Fermi 
energy) consisted of more than 800 basis functions (more than 130 APW/atom) plus Sm 
{5s, 5p), Pd {4p) and Al {2p) local orbitals. The results of our experimental analysis suggest 
that the 4f electrons are localized atomic-like. Therefore the Sm 4f states were treated in 
an open-core approximation with the stable atomic configuration 4f- The wave functions 
in the AS region were expanded up to / = 12 and a density GMAX = 14 was used for 
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the interstitial charge. Brillouin-zone (BZ) integration was performed with the tetrahe- 
dron method^ on a 270 special Appoint mesh (4600A;-points in the full BZ). We carefully 
tested the convergence of the results presented with respect to the parameters mentioned 
and found them to be fully sufficient for all the presented characteristics of the SmPd2Al3 
compound. We also analyzed the temperature evolution of the magnetic entropy. In addi- 
tion we included spin-orbit coupling into the calculations for the electronic structure of the 
valence states (delocalized Bloch states); we found only minor changes in the DOS curves 
and calculated crystal field parameters. The calculations first principles of the crystal field 
interaction were performed using the method described in ref.— . In these calculations the 
electronic structure and corresponding distribution of the ground state charge density was 
obtained using the full-potential APW + lo method. The CF parameters originated from 
the aspherical part of the total single-particle DFT potential in the crystal. To eliminate 
self-interaction, a self-consistent procedure was first converged with the y^/ electrons in the 
core^, which was the open-core approximation used in this work. To diagonalize the mi- 
croscopic CF Hamiltonian of the hexagonal point group symmetry, we used J = 5/2, J = 
7/2 and J = 9/2 multiplets and Wybourne parameterization of the CF matrix elements^. 
Therefore the J-mixing and the intermediate-coupling many-particle ^/-wave functions were 
properly taken into account. The CF perturbation matrix had a dimension of 24 and the 
resulting eigenenergies and eigenfunctions were used to calculate the magnetic susceptibility 
ClS db function of temperature. 



III. RESULTS 



We successfully prepared the SmPd2Al3 single crystal of high quality as confirmed by the 
Laue patterns. The XRPD data contained only the reflections compatible with a hexagonal 
structure of the PrNi2Al3 type(space group P6/mmm) with lattice parameters a = 5.293 A, c 
= 4.064 A, which are in good agreement with published data. The proper 1:2:3 stoichiometry 
of the grown crystal was confirmed by EDX analysis within the accuracy of the method. No 
spurious phase has been located. We also carried out an area analysis and concentrations 
gradients were also studied. No concentration gradients have been detected. 

Specific heat 

The Cp vs. T data were collected by two series of measurements in a magnetic field 
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applied along the a- and c-axis, respectively. In zero magnetic field a large and sharp /-shape 
anomaly in the Cp{T) dependence, peaking at 12.4 K was detected and three much smaller 
peak-like anomalies located at 4.4, 3.9, 3.4 K, respectively (see FigJT]). We tentatively denote 
these temperatures as Tc, Ti , T2 , T^, respectively. The 12.4-K anomaly became gradually 
enhanced with increasing magnetic field applied along the c-axis{B\\c), though it remained 
'pinned' at nearly the same temperature (it move only very slightly to higher temperatures). 
Although such evolution is rather unusual, the magnetization and resistivity data presented 
below confirm that this specific-heat anomaly is apparently associated with the onset of 
a ferromagnetic ordering in SmPd2Al3 and 12.4 K is the Curie temperature (Tc) of this 
material. The low-temperature anomalies, which were smeared out in fields B\\c higher 
than 1 T, were presumably associated with an order-order magnetic phase transition. The 
12.4-K peak remained nearly intact in a magnetic field up to 9 T applied along the a-axis 
(not shown in the figure), thus confirming a strong magnetocrystalline anisotropy in the 
ferromagnetic state. On the other hand, the low temperature group of anomalies became 
continuously smeared out with an increasing field B\\a (not shown in the figure). This result 
indicates that the low temperature phases (T < Ti) have nonzero components (probably 
antiferromagnetic) in the basal plane of the hexagonal structure. For proper comparison, in 
the following discussion of our results along with existing literature data on SmPd2Al3, we 
have collected together all existing information (including our results) on the temperatures of 
anomalies of specific heat, magnetization in very low fields, AC susceptibility and electrical 
resistivity in the Tab. [B It is immediately apparent that there is an unambiguous response 
of all measured properties to the onset of magnetic ordering slightly above 12 K. whereas 
the situation is not so straightforward in the of features observed between 3 and 5 K. The 
temperature dependence of the specific heat of SmPd2Al3, which is presented over a wide 
temperature range in Fig|2l was considered to be a sum of the electronic Cei, phonon Cph, 
magnetic Cmag and Shottky Csch contribution, respectively: 

Cp = Cel + Cph + Cmag + Csch (l) 

The most significant is the phonon contribution mainly at higher temperatures and this 
contribution was evaluated within the afore mentioned Debye and Einstein model^-, the 
results being listed in Tab. [Tll The values of the phonon contribution have then been 
subtracted from the raw specific-heat data. As input data for the used model, the values 
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of 7 = 6.0 mJ/molK^^^, Td = 180 K— and figures of degeneration of the 4 phonon optical 
branches 2-5-3-5^, were taken. The remaining values were analyzed by testing all the free 
parameters in order to reach the best agreement with experimental data, based on the 
equation listed hereafter for the other specific heat contribution: 

1 3N-3 1 

c..-RiY^^c.+ Y.r^^c.) (2) 



Cei = 7T (3) 

The value found for the Sommerfeld coefficient was 7 = 7 mJ/mol K~^, which is a typical 
figure that may be expected and measured in Sm intermetallics compounds^. We also 
analyzed the temperature evolution of the magnetic entropy (Figj2]). In the temperature 
range 1.8 - 20 K we observed a high jump up to the value Rln2 and we presume that the 
magnetic ground state of our compound is a doublet. An additional increase of temperature 
led to a gradual growth of magnetic entropy. The values for entropy should reach the value 
Rln6 when the energy of the third doublet is achieved, but unfortunately this doublet is 
above room temperature. We have measured specific heat data up to 400 K, but the data 
were strongly affected by the transformation of apiezon from temperatures over 250 K and 
it was inappropriate to solve the absolute value of entropy from the experimental data in 
this temperature range. 

A huge energy gap of 130 K was observed between the ground and the first excited 
state, respectively, which is in good agreement with the findings in'^ and our theoretical 
calculations. The existence of a third doublet, indicated to be around room temperature, 
motivated us to measure magnetization up to 400 K, because a strong infiuence of CF on 
the susceptibility behavior was expected, as will be discussed later. 

AC susceptibility and magnetization 

The temperature dependence of the susceptibility measured in the AC magnetic field 
applied along the a- and c-axis, respectively, is shown in FiglHl The a-axis data were 
characterized by a very low, nearly-temperature-independent signal free of any considerable 
anomaly exceeding signal noise, whereas the data measured in the AC field applied along the 
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c-axis exhibited three anomahes, presumably indicating magnetic phase transitions. When 
coohng, the c-axis data exhibited a sudden upturn commencing at Tc = 12.4 K, a sharp 
peak at 3.7 K and a shoulder at 4.5 K. The latter two features may be tentatively associated 
with the magnetic phase transitions at T2 and Ti, respectively, which is tentatively implied 
on comparing the data in Tab. HI 

The theory~, however, says that T^v, the temperature at the maximum of the specific 
heat and the maximum in d{x T)/T vs. T frequently does not coincide with the maximum 
in X vs. T in real systems. That was also the case in our study where the maximum in 
d{xT) / T vs. r was at 3.4 K, which coincided with the Tn value determined from specific- 
heat data. The dramatic difference between the magnetization curves measured at 1.8 K 
in the magnetic field applied along the a- and c-axis, respectively, (FigH]) corroborated, in 
agreement with the specific-heat and AC susceptibility data, the conclusion that the mag- 
netocrystalline anisotropy in SmPd2Al3 was uniaxial and the easy-magnetization direction 
was the crystallographic c-axis. Whereas the magnetic moment in the field B\\c is dominant 
and saturates above 2 T, a very weak paramagnetic-like linear (in higher fields) response 
is observed for B\\a. The slight non-linearity and hysteresis observed in the loop measured 
along the hard magnetization axis was considered to be no more than the result of an imper- 
fect geometry in the experiment and slight mis-orientation of the sample than some intrinsic 
hard axis behavior. The c-axis magnetic moment of 0.19 /i^/f-u. at 1.8 K observed in a field 
of 7 T is nearly Stimes lower than the expected ordered moment for the ground-state multi- 
plet of the free Sm'^"'" ion {gJfis = 0.71 fis), which is presumably responsible for the entire 
magnetic behavior of SmPd2Al3. We roughly estimated the observed magnetic moment of 
the ground state from crystal field calculations. We obtained a semi-quantitative agreement 
with the experimentally-observed magnetization at 1.8 K. For simplicity, we did not include 
the exchange interaction due to its complex nature. The step-like shape of the c-axis curve 
in fields weaker than 2 T was fully reproducible and is apparently connected with the specific 
features of samarium magnetism in SmPd2Al3. A detailed view of the magnetization loops 
for B\\c is presented in FigJSl 

A series of four metamagnetic-like transitions can be identified on the 1.8-K virgin mag- 
netization curve. The onsets of transitions denoted by arrows can be found at 0.03, 0.35, 
0.5 and 0.75 T, respectively. The transitions in lower fields can still be traced on the curve 
measured at 2.5 K, but they have already gone at temperatures 3.4 K and higher. The 
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magnetization loops measured at 1.8 and 2.5 K exhibited not only some metamagnetic-like 
features but also pronounced hysteresis. All these features vanished at temperatures from 
3.4 K. This implies that the ground-state phase, which forms at is antiferromagnetic 
and apparently rather complicated. The magnetization curves measured at < T < Tq, 
however, still showed a strong tendency to saturate as usually found in ferromagnets. The 
Arrott-plot analysis of the magnetization curves confirmed the value of Tc = 12.4 K. Above 
this temperature, the typical paramagnetic response of the magnetization to the magnetic 
field applied parallel to the c-axis was observed. In fields parallel to the a-axis, the magneti- 
zation remained and very weak and nearly temperature independent below Tc irrespective 
of the dramatic developments in the c-axis magnetization. We also measured the tempera- 
ture dependence of the magnetic moment in several constant magnetic fields applied along 
the c-axis, the results of which are shown in Fig|6l In very low fields the value of the mag- 
netic moment sharply increased below Tq, which was in a good agreement with the anomaly 
observed in the specific-heat and AC-susceptibility data as well as with the evolution of the 
magnetization loops with temperature. The 10-mT data were in a reasonable agreement 
with the AC susceptibility behavior. The evolution of 3.4-K anomaly with the increasing 
magnetic field, as well as the evolution of the magnetization curves with temperature de- 
creasing below 5 K, resembled a transition to a low-temperature antiferromagnetic state, 
which is stable only in very low magnetic fields. The bifurcation of the thermomagnetic 
curve in the ZFC and FC branches appeared at temperatures below 3.6 K (Fig|7]), where an 
hysteresis of the magnetization curves also apeared . The bifurcation point shifted to lower 
temperatures with increasing applied magnetic field. However, an alternative interpretation 
can also be considered. It involves a ferromagnetic component to the magnetic order at the 
lowest temperature. If all four zero-field ordering transitions are linked to the field-induced 
transitions, then the lowest-field one may be linked to the ferromagnetic transition at 12.4 
K, and would be in fact the coercive field at which the sample becomes a single domain, 
as far as the ferromagnetic component is concerned. The developing difference (FigJT]) in 
zero-field-cooled and field-cooled magnetization can be considered direct evidence of devel- 
oping coercivity, which in turn is linked to ferromagnetism. The upper three field-induced 
transitions at 1.8 K may then be directly connected to the zero-field transitions at 3.4, 3.9 
and 4.4 K. The internal molecular field can be estimated at 1.8 K as 15 T. This value is 
roughly comparable with the coercive field. A specific issue of magnetism in SmPd2Al3 is 
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the involvement of the populated CF levels of the first and the second exited multiplet. This 
is well documented by the extremely broad Schottky contribution to the specific heat (see 
FigJ2]), and, in addition, a necessary impact on the magnetization (susceptibility) behavior 
at elevated temperatures would be expected. The temperature dependence of the inverse 
susceptibility B/ M is strongly nonlinear for both field directions, i.e. it does not follow the 
Curie- Wiess law. The data in FigJH however, can be well fitted by the modified Curie- Weiss 
law^— between 15 and 150 K: 

C 

X = Xo + (5) 

where the temperature independent term xq ^lay be attributed to the temperature in- 
dependent van Vleck contribution which is due to the small energy difference between the 
ground state multiplet J = 5/2 and the first excited multiplet J = 7/2. The fit yields the 
values for the paramagnetic Curie temperature Op = - 21.3 K and effective Sm moment /ie// 
= 0.87 fiB, which is in a good agreement with the value of 0.85 fiB expected for the Sm^+ 
free ion; Xo = 1.15-10^® m^/mol. The considerable deviation from the fit line at higher 
temperatures can be attributed to population of the second excited multiplet. 

To see the impact of magnetism on electrical transport we measured also the temperature 
dependence of the electrical resistivity for current along c (see Figj9]). We have detected a 
clear cusp around 12.5 K, which is in good agreement to the onset of ferromagnetism at 
Tc determined by other experiments discussed above and also with results presented in the 
literature^'^''. Another cusp-like anomaly is seen around 4 K which coincides with T2. 

First-principles and CF calculations 

The calculated GGA electronic density of states (DOS) of SmPd2Al3, using the exper- 
imentally determined lattice parameters, is shown in Fig JTOT a). The occupied part of the 
DOS has a width of 9.5 eV. The first region, from -9.5 eV to -4.8 eV, consists mainly of 
the free electron-like states from the interstitial region and the Sm-fo, Pd- 5s, Al- 3s and 
Al-3p states from AS spheres (see Fig{TO]b-d). The following band group from -4.8 eV to 
the Fermi level represents mainly the Pd-4d states hybridizing with the Sm.-5d and Al- 3p 
states. The unoccupied states above the Fermi level have predominantly the Sm- 5d char- 
acter with an admixture of the Pd-^d and Al- 3p states and the large contribution from the 
free electron- like interstitial region (FigJTO] b-d). The energy position of the localized Sm 
4f states is correctly below the Fermi level. We performed also the spin polarized LSDA 
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and GGA calculations in order to estimate the value of the hybridization induced Pd and 
Al magnetic moments {4f in core and 4f in band as a Bloch states) and found values less 
than 0.1 fiB and 0.01 fis, respectively. The Fermi level for SmPd2Al3 situated at the de- 
scending part of the DOS yielded N{Ep) = 2.51 states l/eV(f.u.). The orbital analysis of 
the DOS shows that mainly the Pd-^d, Pi-Sd and Al-3p states that contribute to the total 
DOS at Ep- The value of the DOS at Ep is too small to cause a spontaneous magnetic 
polarization of the Pd-4d states. The value of the DOS at Ep for SmPd2Al3 corresponds 
to an electronic specific heat coefficient 7 = 5.9 mJ mol~^ K~^, which is somewhat lower 
than the 7 value of 7.0 mJ mol^^ K^^ derived from our specific heat data. This points to a 
rather low value of the mass-enhancement coefficient A = 0.18 for SmPd2Al3 {-fexp. = (1 + 
A) lband)i indicating a weak electron-phonon interaction in the SmPd2Al3 compound. Next, 
the magnetic susceptibility was calculated on the basis of the CF scheme. The microscopic 
CF Hamiltonian for the Sm atomic configuration 4f in the hexagonal symmetry has four 
independent parameters A^. Reliable CF parameters were obtained for the NdPd2Al3 com- 
pound by fitting inelastic neutron data io^>^; the values obtained were A\ = -386 K, = 
42 K, Aq = 6.8 K and A^ = -134 K. At first we decided to check these values by our first 
principles calculations. 

The first principles calculation of the crystal field parameters leads to the following values: 
Al = -307 K, Al = 25 K, A° = 1.9 K and Al = -88 K for SmPd2Al3 using the GGA forn>i^ of 
the exchange-correlation potential. When we consider spin-orbit interaction the values were: 
Al = -308, Al = 25 Aq = 1.9 and Al= -88. The second-order CF parameter has the correct 
sign, which determines the easy c-axis of SmPd2Al3 at low temperatures in agreement with 
the analysis of our susceptibility and magnetization data. The other parameters have also the 
correct sign and similar values compared to the CF parameters resulting from the analysis of 
inelastic neutron data in^^. Moreover the calculated gap between the ground and the first 
excited doublet is more than 100 K, which is in goof agreement with the Schottky specific 
heat analysis and this fact is also supported by the reliability of both the CF parameter sets 
obtained above. Since our theoretical approach is semi-quantitative only^fl, we may state 
that a consistent description of the CF interaction in SmPd2Al3 compound is obtained. We 
would like to emphasize that we calculated a similar temperature dependence of magnetic 
susceptibility, which has an intersection of curves along the c- and a-axis (see Fig JTTj) . as we 
found by experiment. The intersection was not obtained using the CF acting on the ground 
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state multiplet J = 5/2 only so it requires at least the CF acting on J = 5/2, J = 7/2 and 
J = 9/2 quantum subspace. Note that no intersection has been obtained from the theory 
published by Liu^. 

A detailed comparison would require a correct description of the exchange interaction 
in the paramagnetic region above 12.6 K, since the simple molecular-field and CF model 
was found to provide only a semiquantitative agreement with our experimental data, and is 
therefore definitely too crude for a description of our measured paramagnetic susceptibility 
versus temperature data. Furthermore, we therefore also decided not to fit the measured 
data by a too-crude molecular field and CF model approach. 



IV. CONCLUSIONS 



We have prepared a single crystal of the SmPd2Al3 compound in order to study specific 
features of Sm magnetism. The magnetization, AC susceptibility and specific-heat measure- 
ments of the crystal exposed to magnetic fields applied along the principal crystallographic 
axes revealed a strong uniaxial magnetocrystalline anisotropy (even in paramagnetic state) 
with a easy magnetization direction in the c-axis, which is in contrast to the easy-plane 
anisotropy reported for other compounds in the i?i^Pd2Al3 group. The four magnetic tran- 
sitions observed with temperature dependence of the specific heat at temperatures Ts = 3.4 
K, T2 = 3.9 K, Ti = 4.4 K and Tc = 12.4 K, respectively, and detected in part also in 
magnetization, AC susceptibility and electrical resistivity data, point to a complex magnetic 
phase diagram for SmPd2Al3. Although this material becomes ferromagnetic below Tq = 
12.4 K, an antiferromagnetic ground state seems to become established at low temperatures. 
The series of four metamagnetic transitions detected at 0.03, 0.35, 0.5 and 0.75 T, respec- 
tively, underlines the complexity of the Sm magnetism, which is characterized by a small 
Sm magnetic moment and a complex interplay of crystal-field and exchange interactions. 
The principal role of CF interaction has been confirmed by ab initio calculations and our 
comparison of calculated and experimental paramagnetic susceptibility data. To prove the 
validity of our scenario on the multiphase magnetic diagram, as well as on the magnetic 
phase transitions indicated in our work suitable microscopic experiments are strongly de- 
sired; namely, neutron diffraction and //SR spectroscopy. The high neutron absorption by 
Sm for a standard thermal neutron wavelength implies that a hot neutron beam would be 
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necessary for a successful neutron diffraction experiment with the SmPd2Al3 crystal. 
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TABLE I: Summary of temperatures of anomalies in specific heat, AC susceptibihty, low-field 
magnetization and electrical resistivity observed in our results in comparison with the parameters 
of magnetic transitions in SmPd2Al3 reported in the literature. 



property Source Anomaly location [K] 
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TABLE II: The values and degeneracy of branches of the phonon contribution to the total specific 
heat evaluated on the basis of Einstein and Debye models. 



Branches Degeneracy 9 [K] a [10 ^] 
9d 191 5.0 

Oei 2 103 5.0 

9e2 5 210 4.5 

Oe^. 3 154 4.0 

Oea 5 420 3.0 
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TABLE III: Energy 
J = 5/2 multiplet. 



of the 3 Kramer's doublets after CF splitting of the degenerated ground state 



Levels Energy [K] 
Ai 
A2 130 ± 15 
A3 300 ± 40 
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FIG. 2: (Color online) Comparative analysis of specific heat: the graph (a) shows the total specific 
heat and individual contributions, respectively. Temperature dependence of magnetic entropy (b) 
is shown, as well. 
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FIG. 3: (Color online) Temperature dependence of AC susceptibility measured with the excitation 
field applied along the a- and c- axis. The inset shows detail at low temperature and critical 
temperatures are emphasized by the arrows. 
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FIG. 4: (Color online) Hysteresis loops measured at temperature 1.8 K with applied magnetic field 
parallel and perpendicular on the basal plane. 
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FIG. 6: (Color online) Temperature dependence of magnetic moment of SmPd2Al3 measured in 
selected magnetic fields applied along the c with increasing temperature from the low-temperature 
ZFC (zero-field cooled) state. 
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FIG. 7: (Color online) FC and ZFC thermomagnetic curves measured for SmPd2Al3 measured in 
selected magnetic fields applied along c. 



7 

6 

o 4 



o 



OQ 



3 
2 
1 










= 1 . 1 5e m ^/mol 








^ =-21.3K 

p 


J 


data 




model 











100 



200 
7(K) 



300 



400 



FIG. 8: (Color online) Temperature dependence of inverse susceptibility {B / M) in magnetic field 
applied along the c-axis, respectively. The full line represents the Curie- Weiss fit. 
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FIG. 9: Temperature dependence of relative electrical resistivity of SmPd2Al3 with electric current 
applied along the c-axis in temperature region of magnetic phase transitions. 
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FIG. 10: Total DOS (a) and atom-projected DOS (b-d) of SmPd2Al3. The projected Sm DOS 
(b, thick line), Pd (c, thick line), Al(d, thick hne) and the interstitial region (d, dashed line) are 
shown. Fermi level is put at zero energy. 
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FIG. 11: Paramagnetic susceptibility as a function of temperature of SmPd2Al3 calculated using 
CF model (see text). Calculated curves for the CF parameters from NdPd2Al3 (a) and for our 
first principles calculated CF parameters for SmPd2Al3 (b) and calculation the CF acting on the 
ground state multiplet J = 5/2 only (c). Experimental data are displayed in the bottom panel (d). 
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